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A study of the marker-register dial-tone delay problem in No. 5 crossbar 
switching machines led to a special type of cyclic queuinq model. In this 
paper, we present a method for calculating approximately the steady-state 
delays of an arriving customer. When applied to the marker-register 
problem, the model emphasizes the order in which markers are assigned 
to waiting calls and the fact that part of the markers' time is unproductive 
when an "all registers busy" condition occurs. Some numerical results are 
presented, which agree with the observed phenomenon that, for a constant 
marker load, the delays of Touch-Tone® calls are influenced by the load 
on the dial-pulse originating registers, and vice versa. The results are 
compared to those of a simulation of the same problem. The numerical 
results compare favorably in the range of loads that produce a dial-tone 
speed of between 0.05 and 0.15. 

I. INTRODUCTION 

The queuing model described in this paper resulted from a study of 
the marker-register dial-tone delay problem in No. 5 crossbar switching 
machines. A number of queues with Poisson arrivals of equal rates 
are served in a cyclic order by a server with constant service time. 
Upon arriving at a nonempty queue, the server chooses a customer 
from the queue at random. After one service time, the customer either 
leaves the system with a certain predetermined probability or rejoins 
his queue. In both cases, the server uses a fixed amount of time and 
moves to the next queue ; thus, at most one customer leaves the system 
following each arrival of the server at a queue. 

Related models were treated by Cooper, 1 Cooper and Murray, 2 and 
Eisenberg. 3 In Refs. 1 and 2, the server either empties the queue being 
served or serves all those present at the queue in its arrival epoch. 
The case of two queues with different arrival rates is treated in Ref. 3. 
In these papers, the Laplace-Stieltjes transforms of the waiting time 
distributions were obtained. Attempts to obtain the distributions for 
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similar models by approximate methods were made by Leibowitz 4 and 
Schay, Jr. 6 

The method presented here approximately calculates the steady- 
state delays of an arriving customer. The approximation is carried out 
by modifying the model to achieve a manageable state space. Next, 
the method is applied to the problem of dial-tone delay in the No. 5 
crossbar switching machine. The order in which markers are assigned 
to waiting calls and the fact that part of the markers' time is unpro- 
ductive when an "all registers busy" condition occurs are emphasized 
in the model. The numerical results presented agree with the observed 
phenomenon that, for a constant marker load, the delays of Touch- 
Tone® calls are influenced by the load on the dial-pulse originating 
registers, and vice versa. The results are compared to those of a simu- 
lation of the same problem. 

Several features of the No. 5 crossbar machine, which may have an 
influence on the dial-tone speed, were excluded. Some of these features 
were investigated in subsequent work by H. A. Guess 6 and are de- 
scribed in more detail in Section XVI of this paper. 

II. THE MODEL 

Let N queues Ei, E 2 , • • -, En, N ^ 2, be given. Customer arrivals 
to each queue constitute independent Poisson processes, all having the 
same rate X. The queues are served by a single server in the following 
way: At each point in time, the server is at some queue. Transitions 
in its position occur at discrete time epochs which are equally spaced 
with periods of duration T. At such time epochs, the server moves 
instantaneously to the next nonempty queue in a circular order, chooses 
a waiting customer in that queue at random, and stays with this 
customer until the end of the period. The served customer then leaves 
the queue with probability p*, or remains in the queue with probability 
q* = 1 — p*. If all the queues are empty at a transition point, then 
no change occurs in the position of the server. 

We assume further that T is small with respect to the accuracy 
with which we want to know the delays, and thus all the arrivals can 
be assumed to occur at the transition epochs and the queuing process 
can be considered in discrete time. Thus, if Xi,k is the number of 
arrivals to Ei at the kth time epoch, then all Xi,k, 1 £ i £ N f k £ 0, 
are independent identically distributed random variables, all having 
a Poisson distribution with mean XT. 

III. THE FULL STATE SPACE 

The state of the queuing system at any time epoch is defined as the 
(N + l)tuple (mi, m 2 , • • •, m N , n), where w, is the number of waiting 
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customers at E if and n is the position of the server before the transition. 
We call the set of all such states the full state space. It is then clear 
from the discussion in the previous section that our queuing system, 
with the full state space, is a stationary Markov chain. Thus, in 
principle, one can calculate transient and steady-state probabilities. 
However, the full state space is much too large for practical computa- 
tional purposes. 

Consider, for instance, the dial-tone delay problem where a typical 
number of iV is 15. Then, even if the system is so underloaded that 
each m , can be restricted to be either or 1, we have 2 15 X 15 ~ 500,000 
states. A natural approach, which we follow in the remainder of the 
paper, is to approximate the behavior of our system with systems 
having a smaller size state space. 

IV. THE "BLACK BOX" APPROACH 

Some important facts about the system can be deduced by consider- 
ing only the total number of customers in all the queues s = mi + 
ma + • ■ • + )n N . It is clear that the system with this single state is 
again a stationary Markov chain. It is, in fact, a discrete version of 
an M/D/l queue, with the added feature that a customer who is held 
by the server returns to the queue with probability q*. Alternatively, 
the service time measured in units of T may be considered as having 
a geometrical distribution with mean 1/p*. The traffic intensity of 
the system is the p = XNT/p"; thus, we have Theorem 1. 

Theorem 1 : A necessary and sufficient condition for the nonsaturation 
of the system is \NT < p*. 

Note that, because of the symmetry, a particular queue in the sys- 
tem is saturated if and only if the system as a whole is saturated; 
hence, Theorem 1 provides a saturation condition for all the individual 
queues. 

Let us denote by A* the probability that a Poisson-distributed ran- 
dom variable with mean XNT attains the value *. The system has the 
following transition probabilities. For s > 0, 

Pr (s -> s') = A* s .- s+1 p* + A* s .- S q* 
and 

Pr(O^s') = A... 
The equations for the steady-state probabilities P„ are then 

P s - = A*,P + E 1 (A*,-, +1 p* + A**-a*)P, s' = 0,1, ■■■. (1) 



1=1 



Equations (1) can be solved recursively, starting with P = 1 — p. 
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One can also calculate the generating function 
P(u) = £ P a W 

8=0 

(p*-g)(l-u) 
uq* + p* — ue a{l ~ u) ' 

where a = XNT. By evaluating P'(l), the expected value of the total 
number of customers in the system s is 

8 ~ 2(p* - a) (2) 

V. THE MODIFIED MODEL 

We now consider our original system with a new state space. A 
state will now consist of a triplet (m, M, n), where m = mi, M = m 2 
-f-w?3 + ' ' ' + ™n, and n is the same as previously, namely the posi- 
tion of the server before a transition. It is clear that the new state 
space is much smaller than the full state space ; however, the Markovian 
property is lost. This can be seen by the following argument: If M 
was positive at time k — 1, and if in the transition between k — \tok 
the server skipped a large number of queues, then those M customers 
were concentrated in the remaining queues; thus, it is probable that, 
in the k to k + 1 transition, a small number of queues will be skipped. 

At this point, we modify our model to make it a Markov chain with 
respect to the new state space. To do this, we need to define one-step 
transition probabilities so that the behavior of the modified model 
will approximate that of the original model. Let the position of a 
customer be the queue number where he waits. Our key assumption 
concerns the probability distribution of the positions of the M cus- 
tomers, given the state (m, M, n). 

For the remainder of this section, let us enumerate the queues 
E2, • • • , En by starting with the queue following the position of the 
server and observing the cyclic order, skipping E\. Next, we make the 
following assumptions : Let M > ; then 

(i) The positions of the M customers are independent, identically 

distributed, random variables. 
(ii) The probability that any one of the M customers will be in the 
ith queue (in the new order) is in{M), where 

ri (M) = KM) + <* ~ j - W*> i=l,2,..,tf-l, 

where b(M) = 1/{N - 1) - [R(M)/2M~\ and determination 
of R (M) is described below. 
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The rationale behind assumption (ii) is that a customer is less likely 
to be in a queue that has just recently been visited. The average differ- 
ence between the number of customers in the first and last queues is, 
according to assumption (ii), M[vi(M) — ir N -i(M)^\ = R{M), which 
should be approximately equal to the expected number of arrivals 
during one full cycle of the server. Thus, 

R(M)~\TR (M), 

where R (M) is the expected number of nonempty queues, given M. 
Finally, we approximate R (M) by 

*c*> = w - 1) [i - (f^t)*] . 

where the right-hand side is the expected number of nonempty queues 
among E 2 , • • ■ , En, if the M customers are uniformly distributed. 
We conclude this section by using the new assumptions to calculate 
some probabilities that will be needed later. 

Let J denote the number of successive empty queues following the 
position of the server (when Ei is disregarded), J = 0, 1, ■ ■ ■ , N — 1. 

The distribution of J depends on M . Let 

Qi(M) = Pr (J ^ j - 1) 
q s (M) = Pr (J = j - 1) 
3= 1,2, ■■-,#. 



Using assumptions (i) and (ii), we have 



TAT-l ~\M 

QAM) = I E tt.wJ 

f if M > 
Qn(M) = i 

1 1 if M = 0, 



j= 1, .... N-l 



and 



Qi(M) = Qj(M) - Q j+1 (M) j = l, • ■ ., N - 1 

q N (M) = Qn(M). 



VI. TRANSITION PROBABILITIES FOR THE MODIFIED SYSTEM 

Given a state (m, M, n), the transition probabilities of the position 
of the server can be expressed in terms of the Q,(M) and qj(M)'s as 
follows. 



CYCLIC QUEUE DELAYS 1737 



For m > 0, 

("?„/_„ (M) if ri > n 

Pr (n -> ri) = i Q„- n+1 (M) if ri = 1 

[o if Kn'g n. 

For m = 0, M > 0, 

fq n .-„(M) if n' > n 

Pr (to -► to') = \ if re' = 1 (3) 

[g^r+n'-n-i if 1 < ri ^ rc, 

and, for ??* = 0, M = 0, 

1*1 if ri - re 
Pr (n -» n') = j 

[o if n' ^ n. 

Let us denote by A, and a, the probability that Poisson-distributed 
random variables with means \(N — 1)T and XT, respectively, attain 
the value i. We can now write the state transition probabilities. 

For m > 0, M > : 

Pr [(m, M, n) -» (m', M', re')] 

ra m '_ m g„'_„(M)[AAf'-M+iP* + Am'-m9*] if ri > n 

= iAM'-MQN+l-n(M)£am'- m +lP* + «»,'-.»?*] if ft' = 1, 

[u otherwise. 

For m = 0, M > : 

Pr [(0, M, n) -» (in', M', ri)'] = a m .lA M --M+iV* + Aw-mQ*] 

\q n ._ n {M) if re' > n 

x\0 if n' = l 

[q ir+n ,-. n -. l (M) if Kn'ln. 

For m > 0, Af - : 

Pr [(;/(, 0, n) -> (m', M', re')] 

fAA/'Ca m /_ m+ ip* + a m »_ m g*] if n' = 1, 

(0 if w' 5^ 1, 

and finally, for in = M = 0, 

fa m ' A m> if re' = n 

Pr[(0,0 1 n)-( W ',M',n')]- (4) 

[u if nV re. 

VII. STEADY-STATE EQUATIONS FOR THE MODIFIED SYSTEM 

If we consider the total number of customers in the system s = m 
-\-M, then it is clear that all the results of Section IV remain valid 
for the modified system. In particular, \NT < p* is the necessary 
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and sufficient condition for the existence of a steady state. In the 
following discussion, we assume that this condition is satisfied. Let 
P{m, M, n) be the steady-state probability of the state (m, M, n) 
satisfying the following equations. 

For n' > 1, m' ^ 0, M' ^ 0: 

P(m', M', ri) 

n'—l m' M' + l 

= L E L am'-mqn'-n(M)£A M '-M+lP* + A M '-M<1*1 

n=»l m=0 Af—l 

XP(m,M,n)+ £ "I?* a m .ZA M '-u+iP* + A M .- M q*l 

n = n' M=l 

Xq N+n '-n-i(M)P(0, M, n) + a M .A M >P(0, 0, n'). (5) 
For n' = 1, m' ^ 0, M' ^ 0: 

A' m'+l W 

P(m',M',l) = Z Z Z AM>-AfQ N+ i- n (M)la m ,_ m+l p* + a m .- m q*l 

n=l m-1 M-0 

XP(m, M, n) + a m »AwP(0, 0, 1). 

VIII. NUMERICAL SOLUTION OF THE STEADY-STATE EQUATIONS 

Equations (5) can be solved either as a system of linear equations 
with the auxiliary equation 

Z P(m, M, n) = 1 

m.M.n 

or by starting with any initial distribution and iterating it through (4) 
until a desired degree of convergence is obtained. 

It seems adequate to adopt the second method since the steady state 
of the modified system is of interest to us only as an approximation 
to the steady state of the original system. Thus, an extensive compu- 
tational effort to obtain an accurate solution to (4) is not warranted. 
A good initial distribution to start the iterations can be obtained in 
the following way. We have 

Z Z P(m, M, n) = P. s = 0, 1, ■ • -, 

m+Af—s n=l 

where the P s were computed by the method outlined in Section IV. 
If we divide the s customers uniformly among the queues and make 
the position of the server random, we get the following distribution : 

Po(m, s — m,n) = -=B lm, - , sj P„ 

s = 0, 1, • • ■ m = 0, ■ ■ • , s n = 1, ■ • ■, N, 
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where B(m, 1/N, s) is the binomial probability of m successes out of 
s trials, with probability of success = 1/N. 

IX. EVALUATION OF THE DELAYS FOR THE MODIFIED SYSTEM 

Let /(???, M, ft, A;) be the probability that a customer arriving at 
Ei will wait k service periods before leaving the system, given that 
on his arrival the system went into the state (w?, M, n). The customer 
stays for at least one service period, so k ^ 1. Notice also that ft? ^ 1. 

Theorem 2: The delay probabilities satisfy the following recursive 
equations: 



and 



where 



f(m,M,n,l) = ^Q N - n+1 (M) 



f(m, M , n, k + 1) = <n + <r 2 

??? ^ 1, M ^ 0, 1 ^ n ^ N; k £ 1, 



ffi = 



ZEE a m >- m qn>-n(M)[A M '-M+iP* + A M '-Mq*~\ 

n' = n+l M'=M— 1 m' = m 

Xf(m', M', n', k) if M > 0, n < N 
if M = or n = N 

°2 = E E A M '-iuQN+i-n(M) a m >- m+ ip* + a m >- m q* \ 

M' = M m' = m-l \_ 7U J 

Xf(m',M',l,k). 

Proof: f(m, M, ft, 1) equals the probability that the server moves to 
E\, that the particular customer is selected for service, and that he 
leaves the system after the service period. Hence, the formula for 
k = 1. 

For k ^ 1, we have 

/(w?, M,n,k + \) = E Pr [(m, ilf, n) -> (???', M', n') 

m'.Af' ,n' 

O the customer stays in E{\f(m', M', n', k). 

Using eq. (4), we get that <xi is the part of the right-hand side cor- 
responding to ft' t* 1 and <t 2 is the part of the right-hand side corre- 
sponding to ft' = 1. 

Theorem 2 provides a method for calculating the delays conditional 
on the state. Let f*(k) be the probability that a customer arriving at 
Ei will wait k service periods before leaving the system, given that 
before his arrival the system was in the steady state. 
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Theorem 3: 

f*(k) = E E E /« M', n', k)P*(m', M', »') k = 1, 2, • • -, 

m ' = l W — n' = l 

where 

1 r n'-i ^'+1 

P*« M', 7i') = r -i— P(m', M', n') - a L E <?»'- n (M) 
I — ao L n-l Af-l 

• (Aijf-Af+ip* + ^Af-Mg*)P(m', M, n) 

and 

1 r n M' 

P*(m',M',l) = - \P(m', M', 1) - a E E Q* +1 - n (M) 

1 ~ lO L n=l M=0 

•[p*P(m' + 1, M, n) + g*P« M, n)]l 

/or m' fe 1, Jf' £ 0, 2 £ n' ^ JV. 

Proof: The theorem is valid if it is shown that P*(m', M', n') is the 
probability of the state (?»', M', n') at the point of arrival of the 
customer, say, point u, for all m' ^ 1, M ' ^ 0, 1 ^ »' S N. The 
probabihties of states with m' = is zero, since there is at least one 
customer in E\. We know that at u — 1 the system was in a steady 
state. The transition probabilities between u — 1 and u, conditional 
that at least one customer arrives at E\ t are obtained from (4) by re- 
placing a, with a*, where 

a, i = 1, 2, • • • 



a, = i 



1 — a 



z = 0. 



The expressions for P*(m', M', n') can now be calculated by operating 
on the steady-state probabilities with the modified transition prob- 
abilities and using the fact that the steady-state probabilities satisfy 
eqs. (5). 

X. VALIDITY OF THE MODIFIED SYSTEM 

The difference between the original and the modified systems is in 
the rules of the server movement. In the modified system, the server 
does not follow the cyclic order. However, to calculate the delay dis- 
tribution, we are interested only in the pattern of the time points 
when the server is in Ex, and that, hopefully, is similar to the cor- 
responding pattern in the original system. The degree of similarity 
is difficult to check, except by simulating the original system. Verifica- 
tion of the "reasonability" of the modified system may be made by 
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checking whether each E n receives the same number of visits by the 
server, and if E\ gets the right amount of expected number of waiting 
customers, i.e., 

E E P{m,M,n)~h n =\,---,N 

m =0 M-0 -tV 

and 

CO CO JV £ 

E m E E P(m,M,n)'VT; J 

m = Af=On=l JV 

where s is given by eq. (2). 

XI. THE MARKER-REGISTER SYSTEM 

In the rest of the paper, we apply the model to the dial-tone delay 
problem in the No. 5 crossbar switching machine. Following are some 
operational features of the No. 5 crossbar switching machines which 
are relevant to the dial-tone delay distribution. 

Calls appear on line link frames (llf) . The dial tone markers (dtm) 
which are not busy are paired to the waiting calls. Under "normal" 
operation, i.e., when several dtms are free, the llfs look for available 
dtms according to a fixed preference order. When all the dtms become 
busy, a gate is closed and the dtms serve first those llfs that contain 
waiting calls at that moment. If an llf has more than one call waiting, 
only one call will be served during the gating period. 

When a dtm becomes idle following the "all markers busy" condi- 
tion, it looks for a waiting call according to the following scheme. Each 
dtm has its own order in which it scans the llfs. Thus, for example, 
when there are four dtms and 60 llfs, the first dtm will scan the llfs 
in the natural order from to 59, the second dtm will start at llf 15, 
go to 59, and then come back to to 14, etc.* 

When a dtm locates an llf with waiting calls, it chooses one of those 
calls and proceeds to look for an originating register (or) for the call. 
The above choice may be considered random for all practical purposes.* 
If the dtm finds a vacant or, then it connects the call to the or, and 
the calling customer gets a dial tone. If no or is available, then the 
dtm releases the call, and it continues to wait in its llf and to bid for 
a dtm. In both cases, the holding time of the dtm is constant and 
approximately equal. We denote this time by T, where T is approxi- 
mately 0.25 second. In fact, this time is approximately 0.21 second in 



* This is the recommended arrangement, although not all No. 5 crossbar entities 
observe it. 

' We omit consideration of the systematic preference for serving calls in vertical 
group 2. 
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the case where no or can be found, but we ignore this difference to 
simplify our model. 

We also assume that the distribution of holding times of the ors 
is negative exponential for a conservative estimate of the delay dis- 
tribution. The arrival of calls to each llf is assumed to be Poisson, 
with the rate being equal for all llfs. We denote the rate for a single 
llf by X. Finally, in many cases there are two types of calls, dial-pulse 
and Touch-Tone, where both types are served by the same dtms but 
require different ors. When there are two types of calls, the ratio be- 
tween their arrival rates is assumed to be the same in every llf. 

XII. A QUEUING MODEL FOR ONE TYPE OF CALL 

The system described in the previous section is quite complicated, 
and it appears that, to model such a system and be able to derive 
numerical results from the model, some simplifying assumptions are 
inevitable. One such model was proposed by W. S. Hayward. 7 Its basic 
assumption is that, in order to be served, a call must find both a 
marker and a register idle. Once the marker and register start process- 
ing a call, they act independently of each other, each having exponenti- 
ally distributed holding times. To solve the resulting state equations, 
Hayward introduced a system with one type of server, which approxi- 
mates the behavior of his model. 

The present queuing model emphasizes the order in which the 
markers are assigned to waiting calls, and takes into account the 
fact that the time a marker spends serving a call is nearly the same, 
whether or not it found a free register. 

First we assume that each dtm serves only those llfs which are of 
high priority on its list. Thus, in the example of the previous section, 
the first dtm will serve only the first 15 llfs, the second dtm will serve 
only the next 15 llfs, etc. Such an assumption is justified under heavy 
load conditions. We denote by N the number of llfs which are served 
by one dtm. 

Next, we assume that each dtm serves its llfs in a cyclic order and 
that, whenever it finds a llf with waiting calls, it serves exactly one 
call. This assumption is asymptotically valid under heavy traffic 
loads, because of the gating procedure described in the previous 
section. 

Finally, we assume that whenever a dtm serves a call there is a 
fixed probability p* that an or will be available and thus that the 
waiting time of the call will end (i.e., the customer gets a dial tone). 
This assumption would hold if the availability of the ors is independent 
of the number of waiting calls, which is clearly not the case. This 
assumption will cause our model to somewhat underestimate the 
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delays, while the first assumption tends to overestimate them. The 
value of p* can be taken approximately as the delay probability of a 
call, given that the arrival process of the calls to the ors is Poisson, and 
therefore it can be computed by the Erlang C formula. 

Thus, we arrive at the model which was described in Section II, 
with the server the dtm. The server chooses a customer (call) from 
the queue at random. After one service time, the customer either leaves 
the system with probability p* or rejoins his queue (that is, waits in 
his llf). The server then moves to the next llf having a waiting call. 

XIII. SOME THEORETICAL AND NUMERICAL RESULTS 

It was shown in Section IV that the occupancy of the dtm is 

\NT 

Hence, a necessary and sufficient condition for nonsaturation of the 
system is \NT < p*. Also, the expected total number of waiting calls 
in the llfs which are served by the dtm is 

5 - m^T) (2 - pp,} - 

Thus, for a fixed occupancy of the dtm, the expected total number of 
waiting calls is a monotone decreasing function of p*. The same is 
true for the expected delay, W, since by Little's formula 

\TN P p* ' 
and so 

A standard measure for the quality of service is the dial-tone speed 
(dts), which is the probability that the call will have to wait three 
seconds or more for a dial tone. Figure 1 presents some computed 
values of the dts for various values of p* and X with p held constant 
at three different values. It is seen that dts is also a monotone de- 
creasing function of p*, for a fixed dtm occupancy. 

XIV. A MODEL FOR TWO TYPES OF CALLS 

Consider now the case of a system having dial-pulse and Touch-Tone 
calls; this is the usual situation in No. 5 crossbar offices today. Let 
the arrival rates from each llf be Xi and Xo, and let the probabilities 
of finding available registers be p* and pi for dial-pulse and Touch-Tone 
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NO. OF DTMs=4 
NO. OF LLFs=60 



DTM HOLDING TIME = 0.25 SECOND 




0.25 0.50 0.75 

p' (AVAILABILITY OF REGISTERS) 



Fig. 1 — dts as a function of p* and X for constant dtm occupancy p (queuing model). 

calls, respectively. The loads on the dtm due to the two types of calls 
are 



Pi = — — and 



P2 = 



\ 2 NT 



Pi 7>2 

To approximate the delays of this system, we consider a system with 
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one type of call, for which X = Xi + X 2 and p = pi -\- pi. The appro- 
priate p* of this system satisfies 



* Xi + X2 

p* = 



Xi/pJ + M/pl 

Now the state probabilities can be computed by the approximate 
method; however, it is necessary to modify the formulas for the 
computation of the delay distributions to obtain those distributions 
conditional on the type of call. 

Let fi(m, M, n, k) be the probability that a call of type % (i = 1, 2) 
arriving at llf No. 1 will wait k service periods before leaving the 
system, given that on its arrival the system went into the state 
(m, M, n). The recursive formulas of Theorem 2 have to be modified to 



_ Pi 

Ji\m, ±rj. f iv, xj 

and 



Mm,M,n,l) =%Q N - n+ i(M) 

fit 



fi{m, M, n, k + 1) = ci + <r 2 
m ^ 1, M ^ 0, 1 ^ n £ N, k ^ 1, 

where <x\ is as in Section IX except that / is replaced by /, and 
A M '-AiQN+i-n(M) — — a m '_ m+ ip* 

a m ,. m (J^l q* + !-;)][/.('»', M>, 1, *)]. 



«r 2 = E E 



+< 

The proof of the validity of the modified formulas is along the same 
lines as the proof of Theorem 2. 

XV. NUMERICAL RESULTS 

Several computer runs were made for a typical large system with 
60 llfs, 4 dtms (N = 15), 100 dial-pulse ors and 50 Touch-Tone ors, 
with both dial-pulse ors and Touch-Tone ors having a mean holding 
time of 13 seconds. T was taken to be 0.25 second in all runs. The 
parameters varied were X = Xi + X 2 , the total input rate per llf, and 
a = X1/X2 (the ratio of the rates of the two types of calls). 

Figures 2 and 3 describe the results for a = 2, that is, when the 
ratio of the loads is the same as that of the ors. Figure 2 describes the 
behavior of the occupancies of the ors and the dtm, while in Fig. 3 
the dts is plotted as a function of X. Figures 4 and 5 present the cor- 
responding results for a = 3. The values of X were chosen to be near 
the point of saturation, i.e., where the occupancy of the dtm ap- 
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Fig. 2 — Occupancies of the dtm and ors as a function of X (queuing model). 



proaches 1. In this range, the dts is sensitive to small perturbations 
in X. 

Figure 6 shows the dependence of the dts on a for a fixed X (X = 
0.158 per second). It can be seen that the quality of service deteriorates 
as a diverges from the neighborhood of the ratio of the number of 
dial-pulse ors /number of Touch-Tone ors (which equals 2 in our case). 
This is consistent with the observation that Touch-Tone delays are 
significantly influenced by the dial-pulse or load for a constant offered 
load to the dtm. 

The results were compared with those of a simulation model which 
we constructed for the system described in Section XI. Tables I and 
II present a comparison between the results of the queuing model and 
the simulation. We compare the intensities of input for which levels 
of the dts are reached between 0.05 and 0.25 for a balanced system 
and between 0.05 and 0.15 for an unbalanced system. Examining those 
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Fig. 3 — dts as functions of X (queuing model). 



Table I — Values of A. for which given levels of DTS are reached, 

for Ai/A 2 = 2 

(a) Dial-Pulse Calls 



DTS 


Queuing Model 


Simulation 


0.05 
0.10 
0.15 
0.20 
0.25 


0.1680 
0.1705 
0.1720 
0.1730 
0.1735 


0.1695 

0.1735 
0.1760 
0.1775 
0.1790 



(b) Touch-Tone Calls 



DTS 


Queuing Model 


Simulation 


0.10 

0.15 
0.20 
0.25 


0.1690 
0.1705 
0.1720 
0.1725 


0.1690 
0.1720 
0.1740 
0.1760 
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Fig. 4 — Occupancies of the dtm and ORs as a function of X (queuing model). 



Table II — Values of A. for which given levels of DTS are reached, 

for Ai/A 2 = 3 



(a) Dial-Pulse Calls 


DTS 


Queuing Model 


Simulation 


0.05 
0.10 
0.15 
0.20 


0.1540 
0.1565 
0.1580 
0.1585 


0.1550 
0.1590 
0.1610 
0.1620 


(b) Touch-Tone Calls 


DT8 


Queuing Model 


Simulation 


0.05 

0.10 
0.15 


0.1575 
0.1590 
0.1600 


0.1600 
0.1625 
0.1640 
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Fig. 5 — dts as functions of X (queuing model). 

tables, we observe that the differences in the corresponding figures 
for the two models in these regions are less than 4 percent of the total 
input. Also, it can be observed that the computed dts grows faster 
in the queuing model than in the simulation. The reason is that the 
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Fig. 6 — dts and dtm occupancy as functions of Xi/X 2 for a fixed X, X = 0.158. 

assumption made in Section XII, that the probability of finding all 
registers busy can be computed by the Erlang C formula, is incorrect 
in this region. Figure 7 compares the dts as computed by the two 
models. Again we conclude that the fit is fair in the "critical" region. 
It would have been desirable to validate the results by performing a 
field trial. Such a trial should consist of measuring the dts for a No. 5 
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Fig. 7 — dts as function of X (queuing model and simulation). 

crossbar machine with high constant loads. However, observing Fig. 8, 
where the cumulative distributions of the hourly dts from the simu- 
lation were plotted, one sees that these distributions have long tails. 
This implies that, to get a reasonably accurate estimate for the dts, 
say, with a standard error of 1 percent, we would have to run the 
trials for around 25 hours while keeping the load constant. This seems 
to be a difficult task. 

XVI. SUMMARY AND CONCLUSIONS 

We presented a method for modifying the original model, as pre- 
sented in Section II, to a model which has a much smaller state space. 
Methods were described for calculating the steady-state distributions 
of the states and of the delays in. the modified model. The model was 
applied to the problem of calculating dial-tone delays in the No. 5 
crossbar switching machine. This was accomplished by making some 
simplifying assumptions about the order of service of the waiting calls 
by the markers. The numerical results were compared to those of a 
simulation, and found to be close on an important range of the dts. 
This gives us a certain amount of confidence that both models are 
valid, which is especially important because of the difficulty in validat- 
ing the models by experimental data, as discussed in Section XV. 
However, the reader should be aware that several features of the No. 
5 crossbar machine, which may have an influence on the dts, were 
excluded. 
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Fig. 8 — Cumulative distributions of the hourly dial-tone speed (simulation). 



A major effect not covered by the present model is the effect of 
horizontal group blocking on dial-tone speed and on dial-tone marker 
waste usage. Recently obtained field data and theoretical studies re- 
ported in a subsequent paper by H. A. Guess 6 have shown that dial- 
tone speed and dial-tone marker occupancy can be appreciably in- 
creased by horizontal group blocking caused by high average line link 
frame loads and also by poor load balance. Consequently, the dial-tone 
speeds associated with a given call origination rate in an actual No. 
5 crossbar office may be higher than would be predicted by our model. 
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